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The so-called moist-convective shallow-water model, which incorporates moist convection 
in a simple albeit self-consistent way is used to analyse how intense localized vortices, 
with distributions of horizontal velocity and relative vorticity close to those observed 
in tropical cyclones (TC), evolve and interact with topography on the beta-plane at low 
latitudes. Instabilities of such TC-like vortices are studied first in the /- plane approxima¬ 
tion, and their development , interplay with beta-gyres and the role they play in vorticity 
redistribution and intensification are then analysed along the vortex trajectories on the 
beta- plane, both in dry and moist-convective environments. Interactions of the vortices 
with an idealized topography in the form of zonal and meridional ridges and islands of 
elliptic form and the role of moist convection in these processes are then investigated, 
revealing rich vortex-dynamics patterns. The results can be helpful in crude analyses and 
predictions of the evolution of the barotropic component of TC, of their trajectories over 
the ocean and during landfall and of related condensation/precipitation patterns. 

Key words: Vortex Dynamics, Tropical Cyclone, Shallow Water Model, Moist Con¬ 
vection, Beta-plane, Vortex- Topography Interaction. 


1. Introduction 

Understanding and predicting the structure and trajectories of tropical cyclones 
(TC) are, obviously, tasks of utmost practical importance. Due to complexity 
of the related phenomena, simplified models are traditionally used to gain an 
understanding of the fundamental dynamical properties of TC. The vertical 
structure of TC is the first to be simplified in such models. Let us mention in this 
context the pioneering papers by Ooyama (1964, 1969) and Charney and Eliassen 
(1964) where two-layer models of the tropical atmosphere were introduced in 
order to describe TC. Such models continue to be used, e.g. (Jones et al. 2009). 
An even more drastic simplification of the vertical structure leads to purely 
two-dimensional (2D), or quasi- 2D shallow water type models, which helped 
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to understand the instabilities and changes of the horizontal structure of TC 
during their lifecycle. For example, 2D Euler equations were used in (Schubert 
et al. 1999; Kossin and Schubert 2001; Menelaou et al. 2012) to explain the 
formation of so-called meso-vortices in the vicinity of the TC cores. Switching to 
shallow-water models (Hendricks et al. 2014; Schubert et al. 2016; Lahaye and 
Zeitlin 2016) allows us to include inertia- gravity waves and related radiative 
instabilities in the picture. Let us emphasize that instabilities of TC could be 
captured already in the adiabatic, “dry” models (Schubert et al. 1999). However, 
as is well known, moist convection is an essential ingredient of the TC dynamics, 
see the review of (Montgomery and Smith 2017). It can be included in shallow 
water models in the simplest way as a sink in the mass equation, e.g. Schubert 
et al. (2016), following the seminal ideas of Matsuno (1966) and Gill (1982b). 
A more systematic, although following the same philosophy method to include 
moisture, condensation, and related latent heat release was applied in the moist- 
convective rotating shallow water (mcRSW) model, respectively, in one- and two- 
layer versions in (Bouchut et al. 2009) and (Lambaerts et al. 2011). The mcRSW 
models were recently applied to the analysis of instabilities of TC -like vortices 
in (Lahaye and Zeitlin 2016; Rostami and Zeitlin 2018). It should be emphasized 
that instabilities of TC are usually studied in the /- plane approximation, for the 
simple reason that in order to conduct a standard stability analysis, one should 
start from a stationary state, which is not the case of monopolar vortices on the 
/3- plane. 

In the following, we continue the long tradition of shallow-water modelling 
of TC by using the mcRSW, which allows us to incorporate, in a minimalistic 
way, the dynamical effects of the moist convection, in particular the related 
enhancement of cyclonic vorticity. We also rely on the fact that including the 
topography in shallow-water models is easy and computationally friendly. We 
thus study the evolution and trajectories of the TC-like vortices over the ocean 
and their landfall in the beta-plane approximation with an idealized topography, 
centring the beta-plane at a low latitude, and using crude parameterizations 
of condensation and evaporation, which are simplified as much as possible. We 
show that, surprisingly, already such a rudimentary version of the model captures 
well many of the observed, or simulated in “big” models, characteristics of the 
behaviour of real-life TC. It also allows us to understand some of their basic 
dynamical properties, and to make rough predictions of trajectories. Needless to 
say that, by construction, the model obtained by vertical averaging is able to 
reproduce, in a coarse way, only the large-scale barotropic features of TC. 

In addition to their tentative utility in the applications mentioned above, 
the results presented below reveal some previously unknown, or not sufficiently 
investigated, facts in the domain of vortex dynamics proper, such as instabilities 
and trajectories of localized vortices of high intensity on the beta-plane, or 
interactions of such vortices with the topography. 

The paper is organized as follows. In Section 2 we present the model, with short 
explanations, and discuss its hurricane-like vortex solutions, which are considered 
as proxies of TC throughout this study. In Section 3 we present the results of 
the linear stability analysis of such vortices on the /-plane, analyse how these 
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instabilities are modified on the /3- plane, and study their influence, along with 
the moist convection, upon trajectories of the vortices. Section 4 is devoted to a 
study of interactions of TC -like vortices with idealized topography and coasts 
(landfall), both in adiabatic and moist-convective environments, and the influence 
of these interactions on their trajectories and intensities. Section 5 contains a 
summary of the results and conclusions. 


2. The moist-convective rotating shallow water model (mcRSW) and 
hurricane-like vortices 

2 . 1 . The mcRSW model, a reminder 

A systematic derivation from first principles of, respectively, one- and two- 
layer versions of the mcRSW model was given in (Bouchut et al. 2009) and 
(Lambaerts et al. 2011). Let us recall that the standard shallow -water model of 
the atmosphere can be obtained by vertical averaging of the hydrostatic primitive 
equations in the pseudo-height pressure coordinates, e.g. (Zeitlin 2018). Moist 
convection can be included in the model by adding a vertically integrated specific 
humidity with a condensation sink, which is then linked to convective fluxes of 
mass and momentum using Lagrangian conservation of the (linearized) equivalent 
potential temperature. For our purposes below the simplest one-layer version 
of the mcRSW model with the addition of a source of moisture representing 
evaporation, as described in (Lahaye and Zeitlin 2016; Rostami and Zeitlin 2017), 
and an incorporated topography will be sufficient. No explicit dissipation of any 
kind is included in the model at this stage, see the explanations below. The 
equations of this version of the model are: 

( <9 t v + v- Vv + /(y)zAv= -gVh, 

{ d t h + V-[v(h-b)] = -'yC, (2.1) 

[ dtQ + ViQv) = -C + E. 

Here x and y are the zonal and meridional coordinates on the tangent plane, V = 
(< d x , d y ), v = (u, v), u and v are the zonal and meridional components of velocity, 
respectively, h is geopotential height (thickness), b(x,y) is the topography, f(y ) 
is the Coriolis parameter, which is constant f = fo in the /-plane approximation, 
and is equal to / = /o + [3y in the /3-plane approximation, and z is the unit 
vertical vector; 0 is specific humidity integrated over the air column, 7 is a 
parameter depending on the underlying stratification, C is the condensation sink 
and E is the surface evaporation source of moisture. They are parametrized as 
follows: 

c= Q ~ Q T (Q-Q S ), E = a-^-{Q s -Q)U{Q s -Q). ( 2 . 2 ) 

T \Vmax\ 

Here, Q s is a saturation threshold, which can be taken to be constant, as was 
the case in the early papers on mcRSW. Its dependence on pressure anomaly 
can be readily included, cf. (Bouchut et al. 2009), which we will do below. In 
principle, the evaporation threshold could be chosen to be different from as 
the evaporation, physically, does not take place through the whole air column. 
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We take it to be the same in the crudest variant of the model, in order to avoid 
the proliferation of free parameters. The formula for the condensation in (2.2) is a 
column-averaged version of the standard Betts-Miller parameterization, which is 
used in general circulation models; r is the relaxation time, which is of the order 
of several hours in the tropical atmosphere. The parameterization of evaporation 
in (2.2) uses the standard bulk formula, cf. e.g. (Katsaros 2001), which is of 
use in hurricane modelling (Cronin and Chavas 2019), where we renormalized 
the transmission coefficient by |v maa ,|, the maximum value of velocity over the 
domain, i.e. the maximum wind in the hurricane-like vortices we consider, to get 
the non-dimensional parameter a regulating the intensity of evaporation. This is 
a free parameter, the only one representing processes in the unresolved boundary 
layer in the present formulation of the model, its value in the simulations will be 
given below; denotes the Heaviside (step-) function, which accounts for the 

switch character of condensation and evaporation. As is natural in the context of 
TC modelling, we consider the vortex motion mostly over the ocean, where 6 = 0, 
and the bulk formula for evaporation (2.2) is well adapted. In the regions where 
6 ^ 0 , which correspond to continental surfaces, the properties of evaporation 
change, and and this parametrization should be changed. We just switch off the 
evaporation in such zones, as a first approach. 

By setting 7 = 0, a = 0, and r —>► 00 which, physically speaking, means 
that the time scale of the dynamical processes in question is much less than r, 
we get the standard adiabatic (“dry” ) rotating shallow water equations, with 
the addition of a passive scalar - the moisture field which decouples from other 
variables and is simply advected by the flow. Notice that in the simplest version 
of the model (2.1) the condensed water vapour just drops off the system. In 
this sense condensation and precipitation are synonymous. Precipitable water, 
precipitation and vaporization can be introduced into the model, following the 
same approach as for the inclusion of water vapour and condensation in ( 2 . 1 ), 
giving an “improved” mcRSW (Rostami and Zeitlin 2018), which we will not be 
using below, again, for simplicity. 

A general analysis, with necessary benchmarks, of the fundamental properties of 
the system ( 2 . 1 ), and of the parent two-layer model, including the conservation 
laws and energy budget, can be found in the above-cited papers. We should 
only stress here that in the absence of evaporation the moist enthalpy m = 
h — b — 7 Q is locally conserved, as can be easily seen by combining the second 
and the third equations in (2.1). The conservation of this quantity replaces the 
mass conservation in the “dry” RSW model and, like h in the latter case, m 
should remain positive, otherwise the phase velocity of the “moist” waves becomes 
imaginary, cf. (Bouchut et al. 2009). This imposes a constraint to be respected in 
numerical simulations. Another important variable for analysis and diagnostics 
quantity is potential vorticity (PV), which is a Lagrangian invariant of the “dry” 
system in the absence of dissipation 

(<9t + v • V)q = 0, q= ^±IM , (2.3) 

where £ = d x v — d y u is relative vorticity. In the presence of condensation, but 
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without evaporation, it is the “moist” PV 


Qm — 


C + f{y) 

h — b — jQ ’ 


(2.4) 


which is a Lagrangian invariant of the “moist” system. In the presence of evap¬ 
oration the moist PV is not conserved, as can be straightforwardly deduced from 
( 2 . 1 ): 

« + v.V)*. = 7j -f t^E. (2.5) 

The term on the right-hand side can be a source or a sink of PV, depending on the 
sign of q m . A dissipation, which is absent in ( 2 . 1 ), but which appears implicitly 
in numerical simulations due to various discretization procedures, also introduces 
a local source, or a local sink of the moist PV, depending on the properties of 
the underlying vorticity field. Numerical simulations with the set of equations 
( 2 . 1 ) that we present below were performed using a second-order high-resolution 
well-balanced finite-volume scheme (Bouchut 2007) in a rectangular domain with 
sponges at the boundaries, which allows us to mostly evacuate fast inertia-gravity 
waves generated by various dynamical processes. No explicit dissipation has 
been added. As was shown in previous works with well-balanced finite-volume 
numerical schemes for shallow-water models, the numerical dissipation, globally, 
mimics the Newtonian one, being concentrated in the zones of high spatial 
gradients. The dependence of the saturation threshold on the pressure anomaly, 
which is proportional to r] = h — Ho, where Hq is the mean value, was taken in the 
form Q s — Qq e~°' 3r] « Qq(1 — 0.3 77) , cf. Bouchut et al. (2009). The resolution and 
the domain size varied depending on the dynamical processes under investigation, 
see below. The values of parameters <a, 7 , and Qq are ad hoc at this level, while 
the value of r is taken to be several time steps of the code, to be consistent with 
its physical value, with the scaling we use below. However, 7 can be absorbed 
into Q, with simultaneous renormalization of r and, hence, can be set to be one 
without loss of generality; a can be also absorbed in Q. With the fixed value 
of 7 we take it to be equal to 0.15, to be consistent with previous simulations 
and benchmarks. We should say that this choice is consistent with the empirical 
value ^ 10 -3 of the dimensional transmission coefficient in the bulk evaporation 
formula (Katsaros 2001), which is also used in hurricane modelling (Cronin and 
Chavas 2019), taking into account that vertical integration and renormalization 
by \v ma x\- In all simulations below the initial Q was uniform, and close to the 
saturation value. 

Let us comment on the property of the system of losing mass whenever the 
condensation is switched on, and the related loss of the “dry” energy = 

f dx dy ^(h — b)^- + 9^-^j (Bouchut et al. 2009). This means that the system 

should be considered as a limit of an infinitely thick upper layer of a two- (or more) 
layer system, with standard mass and energy budgets (Lambaerts et al. 2011). 
As was already mentioned in (Bouchut et al. 2009), an obvious way to restore 
the mass conservation is radiative relaxation, which can be straightforwardly 
introduced in the present model, if necessary, as was done, e.g. in (Rostami et al. 
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2018). We, however, will refrain from doing this for two reasons, again to avoid 
a multiplication of adjustable parameters. The first reason is that de facto in 
the simulations presented below the overall mass loss is less than a tenth of 
a per cent of the initial mass. Notice that the above-mentioned dependence of 
Q s on pressure anomaly prevents an excessive mass loss. The second is that 
the reasonable radiative relaxation time, taking into account that the model, by 
construction, applies to large-scale horizontal motions, would be a global climate 
relaxation scale of the order of couple months, e.g. (Cronin and Emanuel 2013) 
and references therein. Even if we consider the relaxation scale of an individual 
hurricane, which is estimated as 40 — 50 days, cf. e.g. (Chavas and Emanuel 2014), 
the typical time scales of dynamical processes we consider and the duration of 
the simulations presented below are of the order of days, so the influence of such 
relaxation on the results is quite limited , which we checked in sample simulations. 
As should be clear already from ( 2 . 2 ), we purposely do not use any detailed 
parameterization of the boundary layer, although such parameterizations exist 
in shallow-water models (Schecter and Dunkerton 2009). Along the same lines 
as evaporation, a bulk formula for the surface drag can be also straightforwardly 
incorporated in (2.1). We do not introduce it either for the reasons explained 
below. Finally, we should say that there is a way to improve the model by 
allowing for horizontal variations of the mean potential temperature, which gives 
the possibility of representing not only the convective mass flux but also the 
heating due to the latent heat release, as was proposed in (Bouchut et al. 2009) 
and recently realized in (Liu et al. 2020). We leave the use of such an ameliorated 
model for a future work. 

2 .2. Hurricane-like vortices in mcRSW 

We first construct exact stationary-vortex solutions of the mcRSW equations 
in the /-plane approximation. They will be later used as initial conditions for 
simulations on the (3- plane. To this end, we rewrite the equations of the one- 
layer mcRSW on the /-plane in polar coordinates (r, 6): 

\ + (^° + r) * A v + g ^ h = °> 

j d t h+l ( d r (rhu ) + dg(hv )) = - 7 C 

{ d t Q+}: ( d r (rQu ) + dg(Qv)) = -C + E, 

where v = (ur, v 6) is the velocity in polar coordinates, and the Lagrangian 
derivative is ^ + ud r + ^3q. The main observation is that, in the absence 

of evaporation, the axisymmetric azimuthal velocity v{r) and thickness h(r ) in 
cyclo-geostrophic equilibrium 

t ; 2 

-b fov = gd r h , (2.7) 

r 

at zero radial velocity u — 0 , and arbitrary constant or azimuthally symmetric 
Q(r ) ^ Q s give an exact solution of (2.6). In a sense, such a solution is an RSW 
analogue of “dry” tropical cyclones discussed in the literature, cf. (Cronin and 
Chavas 2019) and references therein. In the presence of evaporation, in order 
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to provide an exact solution, Q should stay at the evaporation threshold, i.e. 
Q — Q s , with the choice of the threshold made above. Staying at the evaporation 
threshold corresponds to the immediate relaxation limit, where the condensation 
C = — Q s V-v (Bouchut et al. 2009). With E = £7, and constant Q s the 
axisymmetric divergenceless configurations in this limit are, thus, exact solutions 


of (2.6). 


Coming back to the discussion of extra terms in the mcRSW equations in 
Section 2.1, we see that adding to the right-hand side of the thickness equation 
the global radiative relaxation term of the form — h ~ r H ° with a large, as explained 
above, relaxation time r r , where Hq is the non-perturbed thickness, would lead to 
a slow dissipation of the vortex, and of all of its perturbations, while a relaxation 
to an equilibrium profile verifying (2.7) does not affect the solution itself but 
will slowly damp the perturbations. Adding a bulk surface drug of a standard 
form —K( |v|)v to the right-hand side of the momentum equation would lead 
to a spindown of the vortex solution. In what follows we are interested in the 
long-time influence of the beta-effect and topography on the intense hurricane¬ 
like vortices, so to maintain their intensity, and mimic the real-life TC, we will 
not include these terms in the simulations below. Needless to say that the full 
thermodynamics of real-life TC cannot be captured in mcRSW, so the vortices 
we consider are dynamical proxies of these latter. 

Realistic radial profiles of azimuthal velocity and pressure of TC - like vortices 
can be obtained by piecewise-constant, or smoothed piecewise-constant distri¬ 
butions of relative vorticity (Schubert et al. 1999; Nolan and Montgomery 2002; 
Hendricks et al. 2014; Lahaye and Zeitlin 2016). For many purposes it is, however, 
convenient to have an analytic formula for the azimuthal velocity in terms 
of smooth functions. We use the following expression for the non-dimensional 
azimuthal velocity V (r) of the vortex as a function of non-dimensional radius r 
from the centre, which we will call the “a&c” profile: 


V (r) = (r — ro) a e c ^ r r °^, r > ro, a, 6, c>0 


( 2 . 8 ) 


where a, 6, c and ro are parameters for which we will use the two sets of values 
given in Table 1. The profile (2.8) has been already used in our studies of planetary 
jets and vortices (Rostami et al. 2017, 2018). Corresponding non-dimensional 
radial profiles of azimuthal velocity and relative vorticity £(r) are given in figure 
1 , together with those of non-dimensional thickness anomalies recovered from 
(2.7). The profile in the upper panels of the figure has a wide region of small 
negative vorticity beyond the radius of maximum wind (RMW) while the profile 
at the lower ones has a steeper ascending branch before the RMW and a slower 
descending trend out of it, which is close to velocity profiles in the observations 
of TC (Mallen et al. 2005), but leads to non-zero circulation far from the vortex. 
The results of numerical simulations initialized with these two profiles turn out 
to be similar. In order to have direct control of the amplitude of the velocity, we 
introduce its non-dimensional amplitude e and renormalize V (r) as follows: 


V(r) -> V{r) = e 


V(r) 


max\V(r)\ 


(2.9) 
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Table 1. Parameters of the numerical experiments presented below. ID: identification of the 
experiment. MC: moist-convective environment, D: adiabatic environment. Initial water vapour 
distribution is Q s — 0.01, where Q s — 0.7 is the saturated value. In all experiments ro = 0.1, 
c= 14, a = 0.15. 


\ID 

1 D1 1 

MCI 

| D3 

| MC3 

Dbisl 

| MCbisl 

| Dbis3 

MCbis3 

DbisA 

| MCbisA | 

1 « 

1 3 - 5 1 

3.5 

| 3.5 

1 3 - 5 

| 2.25 

| 2.25 

| 2.25 | 

2.25 

| 2.25 

| 2.25 | 

1 * 

| 0.45 | 

0.45 

| 0.45 

| 0.45 

| 0.25 

| 0.25 

| 0.25 | 

0.25 

| 0.25 

| 0.25 | 

1 6 

1 o.i | 

0T 

| 0.3 

1 °- 3 

| 0.1 

| 0.1 

1 °- 3 1 

0.3 

| 0.4 

1 °- 4 1 

| Ro 

| 3.57 | 

3.57 

| 10.73 

| 10.73 

| 4.8 

| 4.8 

1 145 1 

14.5 

| 19.3 

| 19.3 | 




Figure 1 . Profiles of thickness with indicated Chebyshev points which are used for discretization 
(left panels ), and velocity and relative vorticity (right panels ) of TC-like vortices in RSW. The 
most unstable mode of the velocity profile in the upper row (ID : D3) has azimuthal wavenumber 
2, and that of the lower row (ID : Dbis3 ) has azimuthal wavenumber 3. 


271 As is clear from ( 2 . 8 ) the u abc ” profile is placed at a distance ro from the origin. 

272 It is matched at ro with a linear velocity profile in the interval [0, ro] , in a way 

273 to reproduce the observed approximately constant vorticity cores of TC, and to 

274 have overall continuous velocity and vorticity. 

275 The velocity scale U = 200 m/s and the horizontal scale L = 400 km, which is 

276 typical for TC, are used for non-dimensionalization. The time scale is the vortex 
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turnover time T — L/U ~ 0.5h. As already said, we are working in the tangent 
plane approximation. Being motivated by TC, we place our tangent plane at 
a low latitude, 15°iV, which corresponds to /o = 3.7643 • 10 _ 5 s _1 . The non- 
dimensional /o = Lfo/U ~ 7.5 • 10 -2 , which is, in fact, the inverse global Rossby 
number (Ro) constructed from the scales given above. We, however, will be using 
another definition of the Rossby number: Ro = Fmax/Zo^max, where Rmax is 
the maximum of azimuthal velocity (maximum wind), and L m ax is the radius of 
maximum wind (RMW) which is traditional in the TC literature, e.g. (Lahaye and 
Zeitlin 2016). All of the vortices in the experiments presented in Table 1 have large 
Rossby numbers. We, however will call below those with Ro = (10.37,14.5,19.3) 
’strong’, and the other ones ’weak’, bearing in mind the TC nomenclature. In the 
(quasi-) barotropic mcRSW model, we choose Hq = 8.5 km as the vertical scale. 

With the above-described scaling the non-dimensional equation for the radial 
velocity in the /- plane approximation, which follows after non-dimensionalization 
from ( 2 . 6 ), becomes 

-s4 + '*'4 9 '’' < 2 - io > 

We introduced here a non-dimensional amplitude A of the thickness anomaly 
h = 1 + A 77 ), where Hq is the mean thickness, and the Froude number Fr = 

Uj\JgH(). As in (Lahaye and Zeitlin 2016) we will be assuming in what follows 
that A /Fr 2 ~ 1 . The cyclo-geostrophic equilibrium (2.7) between V(r) and the 
corresponding thickness profile H(r ) in non-dimensional terms becomes 


V 2 

r 


+ /o V = H'(r), 


( 2 . 11 ) 


where the prime here and below denotes an ordinary derivative with respect to 
r. The thickness profile for a given V (r) is thus straightforwardly recovered from 
(2.11) by taking the primitive of the left-hand side. For the “afrc” profile, with 
ro = 0 we obtain 


H{r) = l-f 0 


bv r , 


CL - t - 1 

1 r{ 


O + 1 \ b O+l- 


-)G(cr b , 


2 Zl 

Urmnor* w w 


(2 - 12) 

where = max\V{r)\, G(x,a) = [ 1 /T(a)] t a 1 e t dt and r denotes the 
gamma function. 


3. Instabilities of hurricane-like vortices and their evolution along the 
trajectories on the f3- plane 

We start our investigation of hurricane-like vortices by a linear stability analysis 
in the adiabatic environment of the vortex solutions on the /-plane, which were 
described in Sec. 2 . 2 . We then study how these instabilities are modified on the 
beta-plane and in the presence of moist convection, and how these effects influence 
the vortex trajectories. 









310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 

332 

333 


Hurricane- like vortices in shallow-water model 


11 


0.8 


; 0.6 


6 0.4 
0 


0.2 


-Ch 1=1 
-O- 1=2 

• 1=3 


0.1 0.2 0.3 0.4 

e (Amplitute of the Velocity) 



e (Amplitute of the Velocity) 


Figure 2. Growth rates of the azimuthal modes l = 1, 2, 3 as functions of the vortex intensity. 
Left and right panels correspond to the upper- and lower- row configurations presented in figure 
1, respectively. 


3.1. Instabilities of “dry” TC -like vortices in the f -plane approximation 

Analysis of instabilities of the TC-like vortices below follows the lines of (Lahaye 
and Zeitlin 2016; Rostami and Zeitlin 2017). Taking the “dry” equations (2.6), 
i.e. switching off the condensation (7, linearizing them around the vortex solution 
(2.8), (2.11) and looking for solutions in the form of azimuthal normal modes 

0 u,v,rj ) = ile{[iuo(r),?;o(r),r/o(r)]e* (w_u;t) }, (3.1) 


we get the following eigenproblem for eigenvalues ui and eigenvectors [uo(r),vo(r), r/o(rj]: 



IV/r (/o + 2V/r) - D n 

d r V + V/r + f 0 IV/r l/r 
+ (Fr“ 2 + H )D n Li iv/ r 



(3.2) 


where D n is the operator of differentiation in r, which becomes a differentiation 
matrix after discretization using the Chebyshev collocation points (Wright 1964), 
cf. figure 1. Its numerical solution is obtained with the help of the pseudo-spectral 
collocation method. As the gradient of the vorticity changes sign, cf. figure 1, we 
expect a barotropic instability, according to the standard criteria, i.e. complex 
eigenvalues uo — ojq + zct, where a is the growth rate. By obvious reasons, we 
are mostly interested in the fastest-growing most unstable modes. The results of 
the numerical stability analysis are presented in figure 2, where the growth rates 
of the two most unstable azimuthal modes, as well as that of the mode l — 1, 
which is important in the following, are plotted as functions of the parameter e 
which controls Rmax, cf. (2.9), for the two configurations of figure 1. We should 
emphasize that, as it was already stressed in (Lahaye and Zeitlin 2016), the 
instability is sensitive to the fine details of the vorticity profile. Thus, by increasing 
the slope of the ascending branch of velocity and vorticity in the vicinity of the 
core (lower panel of figure 1) the most unstable mode switches from l = 2 to 
l = 3. The structure of the unstable modes with l = 2 and l = 3 is presented 
in figure 3. More details on the stability analysis of TC-like configurations in 
the framework of mcRSW model, including the dependence of the most unstable 
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Figure 3. Spatial distribution of thickness (colours) and velocity (arrows) (left column) and 
radial dependence of u , v, 77 (right column ), with the critical radii indicated by the vertical lines, 
of the most unstable mode / = 2 of the vortex corresponding to ID : D 3, and of the most 
unstable mode l = 3 corresponding to ID : Dbis3. 

mode on the RMW can be found in (Rostami and Zeitlin 2018). Let us however, 
emphasize, that the results above are obtained in the vertically averaged model 
with the crudest parametrization of condensation and evaporation. Studies of 
instabilities of hurricanes in the framework of full three-dimensional primitive 
equations and in cloud-resolving models, with detailed representation of three- 
dimensional hurricane structure and thermodynamics of the moist air can be 
found in (Schecter and Montgomery 2007) and (Schecter 2018), respectively. 

It is known that the nonlinear evolution of the instability in the /-plane 
approximation leads to homogenization of the vorticity in the centre of TC -like 
vortices, and its intensification, (Schubert et al. 1999; Kossin and Schubert 2001; 
Rozoff et al. 2009; Hendricks et al. 2014; Lahaye and Zeitlin 2016). This process 
corresponds to the evolution and different states of TC cores which were found in 
observations (Mallen et al. 2005). We do not investigate here such evolution on 
the /- plane, as in the above-mentioned studies, and pass directly to the /3-plane 
configuration. 

3.2. Instabilities of TC -like vortices on the ( 3 - plane, and the influence of beta-gyres 

The evolution of intrinsic instabilities of hurricane-like vortices on the /?- plane 
cannot be dissociated from a well and long-time known phenomenon due to 
the /3- effect: development of the so-called beta-gyres (Reznik and Dewar 1994; 
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Sutyrin and Flierl 1994) resulting from the Rossby-wave emission by the vortex. 
Corrections to the process of formation and the characteristics of the beta-gyres, 
initially described in quasi-geostrophic models, were studied in much detail in the 
full RSW model in (Reznik and Grimshaw 2001), although in the small Rossby 
number regime. Due to this phenomenon, the cyclonic vortices start moving in the 
north-west direction in the Northern hemisphere (the so-called beta-drift). The 
influence of the beta-effect upon the trajectories and internal structure of TC-like 
vortices was studied in the barotropic non-divergent model (i.e. 2D Navier-Stokes 
equation on the /3-plane) by Chan and Williams (1987); Fiorino and Elsberry 
(1989), and in the shallow-water model by Willoughby (1994). 

Compared to the previous studies, we have at our disposal a much more 
reliable and higher-resolution simulation tool, and the possibility of including 
the important effects of moist convection. We performed a number of numerical 
experiments on the propagation of TC -like vortices on the /3- plane in the 
Northern hemisphere, and development of instabilities, in the “dry” and moist- 
convective configurations displayed in Table 1. A summary of the results, some 
of which we present below, is that the intense vortices are propelled north¬ 
westward by the beta-gyres, the structure of which is sensitive to the strength 
of vortices, and to the presence, or absence, of moist convection. The meridional 
displacement increases with vortex intensity, and north westward drifting strong 
vortices veer eventually eastward, while at equal intensity the influence of the 
effects of moisture upon the drift is relatively weak. At the same time, the 
asymmetry enhanced by the beta-gyres triggers instabilities which develop in the 
core, leading to formation of secondary “eyewall meso-vortices”, which are well 
documented in the observations of TC, e.g. (Kossin et al. 2002). The formation 
and properties of the meso-vortices are mostly sensitive to the details of the 
vorticity profile and the strength of the parent vortex, and are less sensitive to 
the effects of moisture. During nonlinear evolution, the ring of high PV in the 
vortices with higher Ro breaks down faster than in those with lower Ro. The 
eyewall meso-vortices generated by this process can then wander inside the eye, 
and eventually merge. 

We present in figure 4 a comparison of evolution scenarios of weak and strong 
hurricane-like vortices in the moist-convective environment on the beta-plane 
with ID MCI, MC3, respectively. The vortex profile in both cases corresponds 
to figure 1 (upper panels) with the most unstable mode 1 — 2, although for the 
weaker vortex the growth rates of modes l = 2 and l = 3 are close, cf. figure 
2. In contrast with the studies of the evolution of vortex instabilities on the /- 
plane (Lahaye and Zeitlin 2016; Rostami and Zeitlin 2017), the simulations were 
not initialized with a perturbation corresponding to the most unstable mode, 
the instability having been triggered by the asymmetry introduced by the beta 
-gyres. As follows from the figure, the beta- drift is stronger for the stronger 
vortex (see also below), the instability of the core is stronger as well, and the 
secondary eyewall meso-vortices are more pronounced. Their structure at initial 
stages clearly reflects the dominant l = 2 instability of the stronger vortex, and 
a mixture of l = 2 and l — 3 instabilities of the weaker vortex, in agreement with 
the linear stability analysis of Section 3.1. However, the l = 1 mode, whose growth 
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rate was negligible on the /- plane, cf. figure 2, exhibits a considerable growth 
due to the influence of the bet a-gyres, as follows from figure 5. The mechanism of 
intensification in the core region is clear from figure 4: it is due to the merger of 
meso-vortices. The evolution of weak and strong vortices in the adiabatic (“dry”) 
environment is qualitatively similar to that of figure 4 (not shown). 

An important question is about the role of moisture in the evolution of the in¬ 
tensity of the hurricane-type vortices, cf. (Lahaye and Zeitlin 2016). A comparison 
of the evolution of the peak PV in adiabatic and moist-convective environments 
is presented in figure 6. For the weak vortex, the peak “dry” PV is quasi- 
monotonously decreasing, after a transient growth due to redistribution of PV 
produced by discretization errors, and related readjustment process. However, 
this decrease is attenuated by the effects of moisture. For the strong vortex the 
PV evolution exhibits a period of strong growth associated with the appearance 
of pronounced meso-vortices (the second from the top panel of the right column 
in figure 4), and a second period of growth associated with their merger (the 
third from the top panel of the right column in figure 4), the latter one being 
significantly enhanced by the moist convection. The non-conservation of the peak 
PV is due to numerical dissipation, and to the moist convection, if present, but it 
is rather limited, as the overall change of the peak PV is « 7 % for the weak, and 
^ 10 % for the strong vortex. We thus confirm the result previously obtained in 
the /-plane approximation Lahaye and Zeitlin (2016): the evolution of the vortex 
core instability leads to intensification, and that moist convection substantially 
enhances this process. The modification of the instability pattern on the beta- 
plane due to the growing l — 1 mode apparently does not change the overall 
scenario of intensification. 

The structure of the beta-gyres differs for strong and weak (in the sense 
explained in Sec. 2.2 above) vortices, as follows from figure 7, where their 
evolution is presented in the moist-convective case, the pattern in the “dry” case 
being very close. At the initial stages the beta-gyres arise and slowly intensify. The 
anticyclonic gyre is formed to the north-east and the cyclonic one to the south¬ 
west with respect to the vortex core. The centres of the gyres then start to slowly 
rotate clockwise. Although the individual structure of the beta-gyres is sensitive 
to initial velocity and vorticity profile s, our results on the initial formation 
and evolution of the gyres are in global agreement with previous shallow-water 
modelling of TCs, e.g. (Li and Wang 1994), and also with observational data 
analyses e.g. (Kim et al. 2009). As time goes on, a Rossby-wave tail develops in 
the south-eastern sector with respect to the core of the vortex. Initially, the flow 
between the main vortex and the major gyre is north -westward, but as the gyres 
rotate this flow turns northward. 

Condensation is inherent to the evolution of vortices in the moist-convective 
environment (Lahaye and Zeitlin 2015; Rostami and Zeitlin 2017). While it is 
mostly axisymmetric for TC -like vortices in the /-plane approximation (Lahaye 
and Zeitlin 2015), this is not the case on the (3- plane. We found that the time- 
averaged condensation pattern is clearly asymmetric with respect to the vortex 
centre, as follows from figure 8. At the same time, the water vapour distribution 
roughly follows the beta-gyre structure. As is clear from the figure, the rear-right 
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Figure 4. Comparison of the evolution of weak (ID : MCI, left panel), and strong (ID : MC 3, 
right panel) hurricane-like vortices on the /3-plane in moist-convective environment, as seen in 
the potential vorticity field. Snapshots at t — 40, 55, 80,100[L/C7], from top to bottom. 
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Figure 5. Logarithm of the normalized amplitude of the l = 1 component of the azimuthal 
velocity field calculated in the domain r < 0.5 as a function of time during the evolution of the 
weak ( left panel , ID : D 1, MCI) and strong ( right panel , ID : D 3, MC 3) hurricane-type vortices 
in “dry” (solid) and moist-convective ( dashed) environment. Vorticity profiles correspond to the 
upper panel of figure 1. 




Time [L/U] 


Figure 6. Comparison of the evolution of the maximum of potential vorticity of weak 
(ID : MCI left panel), and strong (ID : MC 3 right panel) hurricane-like vortices on the 
/3-plane in the moist-convective (dashed) versus “dry” (solid) environments. 


445 sector, following the vortex motion, is a zone of the highest condensation, while 

446 the front-left sector is the zone of the lowest condensation. This result strongly 

447 resembles observations of convection in TCs over the western Pacific during the 

448 period from 2005 to 2012 (Yang et al. 2016). 

449 3.3. Trajectories of hurricane-like vortices on the beta-plane and inertia gravity wave 

450 emission 

451 The above-mentioned north-western drift of cyclonic vortices due to the beta- 

452 gyres is robust, and can be established already for the point vortices on the f 3 - 

453 plane (Reznik 1992). In shallow water the drift velocity of low Rossby number vor- 
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Figure 7. Formation of (nonlinear) bet a-gyres during the evolution of the weak ( upper panel , 
ID : MCI), and strong ( lower panel , ID : MC 3) hurricane-type vortices in moist-convective 
environment, as seen in the thickness field beyond the radius of maximum wind. 
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Figure 8. Asymmetric distribution of time-averaged condensation for two hurricane-like 
vortices with different Rossby numbers. ID : MCbisl (left panel) and ID : MCbisS (right 

panel). 


454 tices experiences periodic modulations (Reznik and Grimshaw 2001). As shown 

455 above, we do observe a pronounced beta - drift of the hurricane-type vortices. 

456 Trajectories of vortices of different intensity in adiabatic and moist-convective 

457 environments are displayed in figure 9. The drift velocity in our simulations is 

458 not uniform, as can be inferred from figure 9, where the points are plotted at 

459 equal time intervals - for example, vortices slow down while turning north - but 

460 its variations are not periodic, as in (Reznik and Grimshaw 2001), which is, 

461 apparently, due to the largely different Rossby numbers. 
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Figure 9. Trajectories of hurricane-like vortices with different Rossby numbers in dry (D) and 
moist-convective (MC) environments, as defined in Table 1. Time span is from 0 to 315, at the 
interval 15 [L/U]. 


What is important is that we see a clear-cut difference between the behaviour 
of vortices of moderate and strong intensities. While the former follow the 
trajectories predicted in the literature in the no background flow case, e.g. 
(Adem 1956; McWilliams and Flierl 1979; Chan and Williams 1987; Reznik 1992), 
with the moist convection not substantially affecting the trajectory, and “moist” 
meridional velocity being lower than the “dry” one, the latter behave differently. 
Their trajectories are almost straight for a long time, with meridional velocit¬ 
ies much higher than for moderate -intensity vortices, and “moist” meridional 
velocities higher than the “dry” ones. Remarkably, the trajectories of intense 
vortices turn north, and then north-east at late stages, a behaviour which is 
usually attributed in the observations of TC to the influence of the mean wind, 
although there is none in the present simulations. In order to understand the 
dynamical processes at the origin of these north-east turns we present in figure 
10 the outer structure of an intense vortex at the stage when it starts turning 
northward. As follows from the figure, the main anticyclonic lobe of the associated 
beta-gyre exhibits intensification and significant clockwise displacement. As a 
result, the axis of the dipole which it forms with the main vortex turns north¬ 
eastwards, which explains the change in the direction of propagation. The rate 
of this rotation of the lobe is higher in the “dry” case, which is consistent with 
the behaviour of respective trajectories in figure 9. The evolution of the pressure 
field in the moist-convective environment is significantly affected by condensation 
accompanied by emission of the inertia-gravity waves (IGW), which explains this 
difference. Both condensation and IGW emission are displayed in the two lower 
panels of figure 10. As follows from figure 11, where we diagnosed the associated 
wave activity, the outward-propagating IGW originating in the vortex core are 
amplified with increasing Ro in the moist-convective environment. The above- 
described formation, wandering and merger of meso-vortices during the nonlinear 
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evolution of intense TC - like vortices are at the origin of the peaks of the wave 
activity. 


3.4. Summary and preliminary conclusions 

Before including the orographic effects in the next section, let us briefly sum¬ 
marize the above-described results of the simulations of the lifecycle of hurricane¬ 
like vortices on the /3- plane. The vortex instabilities, which are identified in the 
/-plane approximation, are triggered by the inhomogeneity due to the beta-effect. 
Their development leads to formation of the meso-vortices, which subsequently 
merge and lead to intensification of the vortex at its core. The larger the Rossby 
number of the vortex, the more pronounced are the meso-vortices, and related 
processes. At the same time, according to the known mechanism, beta-gyres are 
developed in the vicinity of the vortex. They lead to the north-western beta-drift, 
which is substantially different for vortices of moderate and strong intensity. In 
the latter case, nonlinear evolution and the related change of orientation of the 
beta-gyres lead to a change of orientation of the trajectory, which eventually 
veers north-eastward. The characteristic vortex-dynamics events described above 
(formation and merger of meso-vortices, change of orientation of the beta-gyres) 
are accompanied by related bursts of inertia-gravity wave emissions by the vortex. 


4. Interactions of hurricane-like vortices with idealized ridges and 
islands 

4.1. Preliminaries 

Changes in structure and intensity of TC travelling over the ocean and en¬ 
countering coasts (landfall) are of obvious high importance, and there exists an 
extensive literature on this subject. On the other hand, interactions of vortices 
with topography is a recurrent subject in the oceanographic literature as well. 
Most of the theoretical advances in our understanding of vortex evolution in 
the presence of topography were achieved using the barotropic non-divergent 
model (called the long-wave approximation in oceanography (Gill 1982a)), which 
in the context of shallow-water modelling, neglects the deformation of the free 
surface. A common approach is based on the asymptotic analysis of the equation 
of conservation of PV in this model at small times, for initially axisymmetric 
vortex configurations. It was initiated by Adem (1956), and then extended and 
developed by many others, e.g. (Smith and Ulrich 1990; Grimshaw et al. 1994; 
Richardson 2000). There also exist numerical simulations using the same model 
with simple topography profiles on the /-plane, e.g. (Grimshaw et al. 1994; 
Richardson 2000; Hinds et al. 2016), and on the beta-plane (vanGeffen and 
Davies 1999), as well as laboratory experiments in rotating tanks (Carnevale 
et al. 1991; Flor and Eames 2002). As can be inferred from these studies, e.g. 
(Grimshaw et al. 1994; Richardson 2000), a monopolar cyclonic vortex in the 
/-plane approximation tends to go first across, and then along the topography 
isolines. At the same time, the vortex perturbation generates topographic waves 
starting to propagate along the topography isolines (Grimshaw et al. 1994). The 
details of the vortex behaviour depend on the parameters of the topography (slope 
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Figure 10. Snapshots of the isobars beyond the radius of maximum wind for the strongest 
vortex (Ro= 19.3) of figure 9 in the “dry” (first row , ID : DbisA) and moist -convective cases 
(.second row , ID : MCbisA). Third row : Spatial distribution of condensation corresponding 
to the second row. Fourth row: Corresponding to the second row convection-coupled inertial 
gravity waves, as seen in the divergence field. Colour bar is saturated in the core region. 
T — 255, 270, 285, 300 L/U. 
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Time Time 

Figure 11. Wave activity measured as the modulus of divergence integrated over the 
region beyond the radius of maximum wind during the evolution of TC -like vortices with 
different Rossby numbers in dry and moist-convective environments. Ro — 14.5 (left panel , 
ID : DbisS, MCbisS), Ro = 19.3 (right panel , ID : Dbis4, MCbisA ), respectively. 

and height), and on the intensity of the vortex. On the beta-plane the situation is 
more complex because of the simultaneous influence of the beta- and topographic 
gyres developing around the vortex, and analysis of PV conservation becomes 
more involved. Yet it helps to interpret the results of numerical simulations, 
e.g. Zehnder (1993) where numerical experiments with a shallow-water model on 
the equatorial beta-plane were performed, together with a qualitative theoretical 
analysis. 

Compared to the existing literature on this subject, we have the advantage 
of being able to use a state-of-the-art high-resolution well-balanced numerical 
scheme, allowing to represent with high fidelity both the effects of vortex stretch¬ 
ing due to deformations of the free surface, and the inertia-gravity and topo¬ 
graphic waves, and also to incorporate the effects of moisture, and to study their 
influence. 

4.2. Interactions of hurricane-like vortices with idealized mountain ridges 

We start our study of the effects of topography by analysing the interaction 
of, first, weak and, second, strong vortices, in the sense of section 2.2 with a 
straight infinite meridional and zonal bumps (idealized mountain ridges). Low 
(bmax = 0.1) and high (b max = 0.2), in units of iLo, bumps of Gaussian shape 
were implemented. 

4.2.1. Weak TC -like vortex 

In figure 12 we present a comparison of the trajectories of a weak TC -like vortex 
encountering a low zonal bump in adiabatic and moist-convective environments, 
and the corresponding evolution of its peak relative vorticity. The recurving of the 
trajectories observed in the figure is correlated with formation of the topographic 
gyres discussed in the above-cited literature. They are clearly seen in the upper 
and middle panels of figure 13, where we present the evolution of relative vorticity, 
and can be anticipated from the conservation of PV. The recurving and rightward 
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Longitude [L] Time[L/U] 


Figure 12. Comparison of the tracks of a weak (ID: Dl, MCI) TC -like vortex (left panel) 
over a Gaussian ridge with bmax = 0.1, and related evolution of the peak relative vorticity (right 
panel) in “dry” (white dots, solid) and moist-convective (black dots, dashed) environments. 
Time span from 0 to 495, the interval between consecutive dots is 15 [L/C/]. 

shift start earlier in the moist-convective environment. This is due to the fact 
that it is the moist enthalpy h — b — 7 Q which replaces the thickness h — b in 
the expression of the conserved (in the absence of dissipation and evaporation) 
potential vorticity, cf. (2.4), which means that the effective topography is higher in 
the moist-convective case. The evolution of the peak vorticity is similar in the two 
cases, although the peak vorticity in the moist-convective environment is higher, 
starting from the encounter of the vortex with the ridge, and this difference is 
increasing in time, cf. the right panel of figure 12. This is due to the condensation 
in the vortex core during the crossing, which is clearly seen in the lower panel of 
figure 13, as condensation enhances the cyclonic vorticity. The minimum of the 
peak vorticity at later times, observed in both curves in the right panel of figure 
12 corresponds to the passage by the TC of the topography maximum, and can be 
easily explained in terms of PV conservation: diminishing the denominator in the 
PV (2.3) should be compensated by diminishing the numerator, i.e. diminishing 
the relative vorticity, and vice versa. A similar form of trajectory, and of the 
peak vorticity evolution, was reported in vanGeffen and Davies (1999), where 
numerical simulations of monopolar vortex interacting with ridges of various 
orientations were performed in the framework of the “dry” non-divergent model. 
Notice that the overall initial depletion of the relative vorticity is due to the 
meridional propagation which increases the planetary vorticity /o + [3y in the 
Northern hemisphere which, again, as follows from the PV conservation, should be 
compensated by a decrease in relative vorticity, at constant thickness. The lag of 
this minimum in the moist-convective case is explained by the faster propagation 
speed in the dry case, which can be inferred from the left panel of figure 12. 
Increasing the height of the ridge leads to a slow-down of the weak vortex over 
the ridge, and to a backward motion in the case of meridional ridge, as follows 
from figure 14. We will come back to this phenomenon in the next subsection, 
and concentrate first on the topographic gyre formation, and its interaction 
with the primary vortex and its beta-gyre, which is illustrated by figure 15. 
The figure displays a rich vortex dynamics pattern with several distinguished 
features. A cyclonic topographic gyre arises by the process of enhancement of 
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Figure 13. Upper row: Spatio-temporal evolution of the relative vorticity during the passage 
of a weak TC -like vortex (ID : MCI) over a zonal Gaussian ridge with bma X =0.1 at y = 3. 
Middle row: Corresponding evolution of h. Lower row: Corresponding evolution of condensation. 
Horizontal dotted lines indicate the ridge. T = 350 corresponds to the minimum of the relative 
vorticity in figure 12. 


the second (cyclonic) lobe of the beta gyre when the TC approaches the ridge. It 
then intensifies. The topographic gyre forms an asymmetric dipole with the rear 
part of the first anticyclonic lobe of the beta gyre, with a tendency of moving 
backward from the ridge along its axis and hence, by conservation of momentum, 
pushing the main vortex in the opposite direction, which can explain the initial 
acceleration of the vortex over the ridge observed in the left panel of figure 14 
(larger distance between the consecutive dots). Recurving of the trajectory is 
due to the deformation and stretching of the anticyclonic vorticity surrounding 
the core, together with that of the near field of the beta-gyre, which leads to 
reorientation of the axis of the vortex dipole which it forms with the main vortex, 
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Figure 14. Trajectories of a weak TC -like vortex (ID : D 1) (white dots indicating peak 
vorticity) passing over a zonal (left panel) and a meridional (right panel) ridges with bmax = 0.2 
in “dry” environment. Time span is [0, 270] and time interval between consecutive dots is 15 L/U. 




and hence a change in the direction of propagation. Stalling of the vortex at the lee 
side of the ridge, which is seen at later stages in the same figure is due to formation 
of a tripolar structure, involving the remnants of the main anticyclonic lobe of 
the beta-gyre, the main vortex, and a secondary anticyclonic topographic gyre, 
which arises downslope and westward of the main vortex. Notice the vorticity 
pattern appearing to the right of the main vortex structure. It is related to the 
topographic wave, which is coupled to the cyclonic topographic gyre, and leads 
to the alignment of the latter along the ridge at the later stages - see also below. 

4.2.2. Strong TC -like vortex 

Figure 16 represents a counterpart of figure 12 for a strong TC -like vortex 
passing over the ridge with b max = 0.2. The trajectories in both “dry” and moist- 
convective environments are looping downslope, leading to trapping of the vortex 
by topography. As follows from figure 17, where we present the corresponding 
evolution of relative vorticity in the “dry” case, a cyclonic topographic gyre 
also arises, and forms a dipole with the tail of the anticyclonic beta-gyre lobe, 
as in the case of the weak TC -like vortex above, with a similar acceleration 
effect, cf. figure 16. The recurving of trajectory at the lee side of topography, 
and subsequent looping are similar to the weak vortex case, and are due to a 
secondary anticyclonic topographic gyre, arising and merging with the remnants 
of the primary beta-gyre and forming a vortex pair with the main vortex and its 
anticyclonic vorticity ring , with a clockwise-rotating axis. Contrary to the weak 
vortex case, no clear tripolar structure is observed, and the cyclonic topographic 
gyre is not aligned with the vorticity field of the topographic wave (see below). 

We should emphasize that in simulations of vanGeffen and Davies (1999) the 
vortex was always crossing the ridge, even a twice higher one , with respect to 
the non-perturbed depth , compared to our set-up, and continued to move in 
the north-west direction, although its trajectory was recurving over the bump. 
Although a direct comparison is not straightforward, as the simulations in vanGef¬ 
fen and Davies (1999) were done with the non-divergent model on the equatorial 
beta-plane with /o = 0, and the vorticity profile and intensity of the vortex were 
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Figure 15. Interaction of a weak TC -like vortex (ID : D 1) in the “dry” environment with a 
zonal Gaussian ridge with 6max = 0.2, as seen in the relative vorticity field (colours). Isolines 
of relative vorticity are zoomed out of the vortex core. Solid red lines: isolines of b. 



Figure 16. Comparison of the tracks of a strong TC -like vortex (ID : T>3, MC3), left panel , 
and evolutions of the related peak relative vorticity (right panel) over the Gaussian ridge with 
bmax = 0.2 in “dry” (white dots, solid line) and moist-convective (black dots, dashed lines) 
environments. Time span is from 0 to 465, and the interval between the dots is 15 [L/C/]. 


630 different, this discrepancy emphasizes the role of divergence effects and related 

631 vortex stretching, which are absent in the non-divergent model. 

632 The evolution of moisture and condensation during the passage of a strong 

633 TC over a zonal ridge is presented in figure 18. The condensation pattern 

634 is similar to that of figure 13, but its intensity at the locations of primary 
ess and secondary vortices is much higher. As follows from the figure, the first 
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Figure 17. Evolution of relative vorticity of a strong TC -like vortex (ID : D 3) passing over a 
zonal Gaussian ridge with frmax = 0.2 (white dots in figure 16). Isolines of relative vorticity are 
zoomed over out of the core areas. 

condensation pattern develops upwind and downslope, but then a second pattern 
moving along the ridge appears. It is related to a topographic wave. The influence 
of the moisture upon the trajectory can be inferred from figure 16. As follows 
from the figure, the trajectories in the “dry” and moist - convective simulations 
are qualitatively similar. As in the previous case of weak TC- like vortex, and 
by the same reason, the trajectory in the moist-convective environment starts 
the rightward shift earlier. The looping is explained similarly, although there is 
a difference of curvatures of ‘dry’ and “moist-convective” trajectories. In view 
of the “dry” evolution, with its stronger curvature, the variations of Q, and 
related condensation (cf. figure 18) apparently contribute to the flattening of the 
trajectory when the vortex exits the ridge. It should be not forgotten that effective 
topography is higher in the moist-convective environment, as already explained. 
The curves of the evolution of peak relative vorticity presented in the right panel 
of figure 16 show a quasi-identical behaviour, apart from the fact that the peak 
vorticity in the “moist-convective” case is higher during the periods when zones of 
condensation are present, as it should be, as condensation enhances the cyclonic 
vorticity. The consecutive minima of the peak vorticity starting at T — 235 
correspond to the repetitive passages by the TC of the topography maximum, 
and, again, can be straightforwardly explained in terms of PV conservation. The 
fact that high topography can trap intense cyclones arriving at incidence angles 
close to 90° is quite remarkable. We will come back to this phenomenon in the 
next section. 

As follows already from the upper panel of figure 18, the passage of TC- 
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Figure 18. Upper row: Evolution of initially uniform moisture during a passage of a strong 
TC -like vortex (ID : MC 3) over a zonal ridge with maximum frmax = 0.2 at y — 3. Lower 
row: Corresponding evolution of the condensation. Horizontal dotted lines indicate the position 
of the ridge. T — 235 corresponds to the first minimum of the relative vorticity in figure 16. 

1 ike vortices over the ridge produces pronounced convergence zones, as it is 
only convergence or divergence of the velocity field which can redistribute the 
initially uniform moisture. The corresponding evolution of the divergence field 
is presented in the upper panel of figure 19, together with its counterpart for 
the meridional bump (lower panel). As follows from the figure, a large-scale 
topographic wave is generated during the encounter of the vortex with the 
ridge, and starts propagating along the ridge. The divergence anomaly due to 
this wave is at the origin of zonally propagating along the ridge patterns of 
humidity anomaly and condensation. As is clear from the previous discussion, 
the topographic wave also influences the dynamics of the vortex system over the 
ridge in the case of the weak TC- like vortex. 

4.3. Interaction of TC-like vortices with an idealized Taiwan-like island 

The Taiwan island, of a roughly elliptic form, with high mountains, is often 
stricken by Western Pacific typhoons having a horizontal scale comparable to the 
scale of the island. This is why observational and theoretical studies of interactions 
of TC with such island are recurrent in the literature, e.g. (Chang et al. 1993; 
Yeh and Elsberry 1993a,b; Kuo et al. 2001; Tang and Chan 2015; Wu et al. 
2015; Tang and Chan 2016; Huang et al. 2016; Lin et al. 2016). The observed 
TC tracks exhibit a variety of patterns with both northward and southward 
deflections. They may be continuous or discontinuous, and can be accompanied 
by formation of secondary lows (Chang et al. 1993; Yeh and Elsberry 1993b; Kuo 
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Figure 19. Isolines of divergence (colours) produced by a strong vortex (ID : MC 3) passing over 
a zonal (upper panel ) and a meridional (lower panel) ridges with bmax — 0.2 in moist-convective 
environment. Small-amplitude signal corresponding to weak inertia-gravity waves produced 
during the encounter, and propagating radially out of the ridge, are filtered out. 


et al. 2001). The process of deflection includes many environmental factors and is 
very complicated, as stressed by Kuo et al. (2001). Many of the modelling studies 
(Yeh and Elsberry 1993a; Wu et al. 2015; Huang et al. 2016; Lin et al. 2016) 
use the simplistic /-plane approximation, with a meridionally oriented island 
and an imposed easterly mean zonal wind which pushes the cyclone towards the 
island. An accent was made on fully developed parameterizations of convection 
and micro- physics, while the hydrodynamical set-up remained over-simplified in 
recent works. The simulations of the TC trajectories near Taiwan on the beta- 
plane with WRF model were performed in (Tang and Chan 2015, 2016). As is 
clear from this short review, most of the modelling in this context was performed 
with simplified models and the seminal study by Kuo et al. (2001) used the “dry” 
shallow-water model. 

In order to understand the role of hydrodynamic processes modified by crudely 
parameterized convection, we undertook a similar study in the mcRSW on the 
beta-plane. We recall that Kuo et al. (2001) demonstrated a possibility of full or 
partial trapping of the cyclone by the island. We can hereby extend such study 
by including together the beta-effect, the topography and the moist convection. 

To represent an idealized Taiwan-like island we implemented a topography in 
the form of a bump of elliptic shape with a Gaussian cross-section, with the 
main axis oriented from south-west to north-east. The elevation of the bump was 
chosen to be smooth, in order to prevent formation of hydraulic jumps and wave 
breaking, with the maximal non-dimensional height b max = 0.2. As is known from 
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Longitude [L] Longitude [L] 

Figure 20. Left panel : Comparison of trajectories of a weak TC- like vortex with 
Ro — 3.57, ID : MCI in the moist-convective environment in the presence (black circles), 
and in the absence (white circles) of the “Taiwan” island with bmax — 0.2 (red isolines). Time 
interval between subsequent circles is 15 [L/U]. Right panel : Same as in the left panel, but for 
a strong TC- like vortex (Ro — 10.37, ID : MC3). 

the above-cited literature, deflection during the approach and (partial) trapping 
by topography are the main features of the TC trajectories near the island. In 
order to check these effects in the mcRSW we initialized the TC-like vortices in the 
vicinity of the island, in a way that their trajectories, which were determined in 
the experiments of section 3.3, cross the location of the island. The initial moisture 
field was uniform, and close to saturation, as before. The numerical experiments 
were performed both with weak and strong TC, at different incidence angles 
with respect to the major axis of the island, and different incidence locations. 
Variations of the incidence angle were achieved by varying the orientation of 
the major axis of the elliptic bump, as we were using pre-calculated trajectories 
for initialization. The results on the track deflection of weak and strong TC- 
like vortices with RMW of the order of the transverse size of the island, and 
trajectories initially approximately perpendicular to this axis are presented in 
figure 20. As follows from the figure, both weak and strong TC- like vortices 
are partially trapped, and start moving along the clockwise path around the 
island, which is in agreement with some of the observed typhoon tracks across 
Taiwan, with the simulations of (Kuo et al. 2001), and also with the results of 
the section 4.2.2 above. The mechanism of recurving of trajectories due to the 
interactions of primary vortex, its beta-gyre, and topographic gyres was described 
in the previous section, and remains qualitatively the same. Its main ingredient 
is the primary cyclonic topographic gyre, which is clearly identified in figure 21, 
where we present the evolution of thickness during the landfall of strong and weak 
TC- like vortices. As is known from the observations, a secondary vortex can be 
formed at the lee side when a TC approaches the island, e.g. (Yeh and Elsberry 
1993a, b) . As follows from figure 21 the secondary low is indeed generated at the 
northern side of the island during the passage of TC. After full intensification 
it detaches but its subsequent evolution depends on the incidence location and 
angle, see below. Such post-landfall phenomenon was documented, for example 
, in the observational study of the typhoon Mindulle (Lee et al. 2008). One 
of the practically important characteristics of the TC landfall on the island is 
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Figure 21. Upper left (right) panel : Selected isobars of a weak TC- like vortex (ID : MCI) 
at T = 210 in the presence (absence) of the “Taiwan” island. Lower row. Same as in the upper 
row, but for the strong TC- like vortex (ID : MC3). 
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Figure 22. Spatial distribution of the accumulated condensation during the passage of strong 
TC - like vortex (ID : MC3, time span from 20 to 270[L/C/], (left panel)), and weak TC - 
like vortex (ID : MCI, time span from 20 to 500[L/C/], (middle panel)) over idealized Taiwan 
island. Right panel : spatially integrated condensation corresponding to the aforementioned cases: 
(ID : MCI, solid, and ID : MC3, dashed). 


732 an amount of precipitation and its evolution in time. Accumulated amount, and 

733 spatio-temporal evolution of condensation produced by the landfall of TC-like 

734 vortices in the middle of the island, at incidence angle ~ 90° are presented in 

735 figure 22. It should be not forgotten that condensation and precipitation are 

736 synonymous in the simplest version of the model we are using. 

737 A panorama of the evolution of all essential fields during an encounter of 

738 a TC - like vortex with the northern part of the idealized Taiwan island at 

739 a smaller incidence angle is presented in figure 23. Typical convergence and 

740 related condensation patterns are seen at the south-western and north-eastern 
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Figure 23. Spatio-temporal evolution of divergence, condensation, relative vorticity, and 
selected isobars of a weak TC - like vortex (ID : MCI) passing over the northern part of 
idealized Taiwan island (red contours), respectively, from top to bottom. 


741 parts of the island. Deformation of the beta-gyres because of the topographic 

742 influence is clearly seen, as well as formation of a secondary low at the northern 

743 tip of the island. The later stages of this landfall are presented in figure 24, 

744 which continues figure 23, and exhibit s a rather complicated vortex system 

745 comprising the primary vortex, the main anticyclonic beta-gyre, which intensifies 

746 and becomes a rather intense localized anticyclonic vortex, and the topographic 
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gyre, which gives rise to a secondary cyclone. Overall, the vortex system evolution 
scenario is similar to the one observed for the straight infinite ridge, cf. figure 15, 
although the terminating topography introduces important new features. Thus, 
the topographic gyre detaches from the island, and together with the primary 
vortex and its anticyclonic beta-gyre forms a three-vortex system moving out of 
the island. At the same time, a secondary cyclonic topographic gyre develops and 
intensifies at the lee side of the island, as follows from figure 24. We should stress 
that the peak vorticity of the main vortex is smaller than that of the detaching 
topographic gyre during some period of time. A documented in the literature 
process of a secondary low developing and replacing the original centre for TC 
(Yeh and Elsberry 1993b; Lin et al. 1999; Wu 2001), as well as disruption of 
trajectories of the observed TC during the landfall, e.g. Kuo et al. (2001), could 
be explained by a process of this kind, as already in our simple model we see a 
swap of the maximum vorticity between the primary and the secondary cyclones. 
The corresponding trajectory of the primary vortex, the related evolution of its 
peak vorticity, together with the peak anticyclonic vorticity, and comparison of 
the peak vorticities of the primary and secondary cyclones are displayed in figure 
25. As follows from the comparison of the left panel of this figure, and figure 
20, trajectories of TC - like vortices are very sensitive to the landfall location 
and incidence angle. The evolution of the anticyclonic vorticity in the middle 
panel of figure 25 clearly shows the emergence of the secondary anticyclone, and 
subsequent re-intensification of the primary cyclone, while the right panel displays 
a swap of the maximum of vorticity between the primary and the secondary 
cyclones. It is worth noticing that our results in the case described in figures 23 - 
25 are similar to those obtained with the far more exhaustive WRF model, and 
also with observational studies (Yeh et al. 2012; Tang and Chan 2016). 

4.4. Resume of the results on TC- like vortices interacting with topography 

The main feature of trajectories of TC- like vortices over the meridional and 
zonal ridges is their recurving, which can lead to looping and trapping. This 
behaviour is due to a rather complex vortex dynamics involving the primary 
vortex, its beta-gyre, which is deformed and stretched during the passage over 
the ridge, and of the primary cyclonic and secondary anticyclonic topographic 
gyres, with an influence of the topographic wave generated during this process. 
The condensation induced by the passage of TC-like vortices over the ridges 
enhances the cyclonic vorticity and displays a specific pattern, which is correlated 
with topographic waves generated during the interaction of the vortex with 
topography, and propagating along the ridge. 

The above-described features also manifest themselves when topography rep¬ 
resents an elongated island. Although the recurving and partial trapping are also 
observed, the terminating topography leads to appearance of “free” secondary 
vortices, arising from the beta-gyres interacting with topography. These vortices 
can have higher peak vorticity that the original vortex, which can explain seem¬ 
ingly paradoxical disruptions of TC tracks over Taiwan reported in the literature. 
The scenarios of the landfall strongly depend on the angle of incidence and its 
location at the island. We should add that many of the detected in this part of the 
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Figure 24. Late stages of the evolution of relative vorticity ( upper row) and PV ( lower row) 
of the landfall of weak TC - like vortex (ID : MCI) continuing figure 23 



Longitude [L] 




Time [L/U] Time [L/U] 


Figure 25. Left panel: trajectory of a TC - like vortex hitting the northern part of the “Taiwan” 
island with & macc = 0.2 (red contours) in the moist-convective environment (black circles) 
compared to that without topography . Time interval between the circles is 15 [L/U]. Middle 
panel: Corresponding evolution of the anticyclonic ( dashed) and cyclonic (solid) peak relative 
vorticities. Right panel : Swap of the maximum of vorticity between the primary vortex (solid), 
and the secondary (dashed) terrain-induced cyclone. 


792 study features are also observed in “dry” simulations with the same topography, 

793 which we do not present. Yet, the moist convection boosts secondary cyclones, 

794 and makes their influence more pronounced. 
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795 5 . Summary and discussion 

796 We have, thus, shown that the mcRSW model produces evolution and tra- 

797 jectories of TC- like vortices, which are qualitatively and quantitatively close, 

798 within the limits of the model, to the observations of TC and the results of 

799 studies with much more sophisticated models. The main advantage of the model 

8 00 is simple and computationally friendly incorporation of crudely parametrized 
soi moist convection, and of arbitrary topography. The simplicity of the model, and 

802 the robustness of modern finite-volume numerical schemes, allows for long-time 

803 high-resolution simulations, and hence for efficient tracking of the intense cyclonic 

804 vortices in different environments. We thus, have shown that TC-1 ike vortices 
80 s in the Northern hemisphere experience the north-western beta-drift, which is 

806 significantly influenced by the moist-convective phenomena in the case of strong 

807 vortices. Our main discovery in this direction is that, in the absence of the 
80 s mean wind, the strong cyclones are travelling much farther North than their 
809 weak counterparts, and experience, after some time, a north-eastern turn, with 
sio a northward displacement being more pronounced in the presence of the moist 
8n convection. Such turns are reported in the observations of TC, but are usually 
812 attributed to the influence of the mean wind. We have also shown that the 
sis intrinsic instabilities of the intense hurricane-type vortices, which were previously 
8i4 identified in the /- plane approximation, are triggered by the beta- effect, and 
sis evolve on the beta-plane as well, although with asymmetric distortions due to the 

8 1 6 presence of the beta-gyres. They lead to formation of secondary eyewall meso- 

8 1 7 vortices which, in turn, evolve, merge, and lead to intensification of the vorticity in 
sis the vortex core. Meso-vortices are recurrently reported in the observations of TC. 

8 1 9 The moist convection enhances the effect of intensification, while the instability 

820 mode with azimuthal wavenumber one, which appears due to the beta-effect, does 

821 not qualitatively change the evolution of the main instability and the dynamics of 

822 meso-vortices. We have demonstrated that the evolution of TC- like vortices leads 

823 to enhanced inertia gravity wave activity in the moist-convective environment, the 

824 peaks in the wave activity being correlated with different stages of the evolution of 

825 the meso-vortices. Concerning the condensation, we found a distinct asymmetry 

826 in the aggregate condensation between north-western and south-eastern sectors 

827 with respect to the vortex centre, again, in accordance with observations of TC. 

828 Using the advantage of easy and numerically friendly incorporation of topo- 

829 graphy in the shallow-water model, we analysed the passage of TC- like vortices 

830 over a ridge-like and mountain island like topographies and demonstrated trap- 

831 ping of the cyclones by topography, accompanied by generation of topographic 

832 waves, which influence the condensation/precipitation patterns. The character 

833 of the trapping (full or partial) depends both on the intensity of the cyclones 

834 and on the characteristics of topography. Our simulations reveal rich patterns of 

835 vortex dynamics, involving the primary vortex, the main anticyclonic lobe of its 

836 beta-gyre, and the topographic cyclonic gyres, which determine the trajectory 

837 and the evolution of the intense vortices. In the context of deflection of the TC 
sss tracks over the mountain Taiwan-like island, we observed several scenarios, from 

839 deflection to trapping, depending on the strength of the cyclone, incidence angle 

840 and landfall location. We found, again, complex scenarios of vortex dynamics 
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of the primary cyclone, a secondary cyclone arising from the topographic gyre, 
and an anticyclone produced by the beta-gyre, with swaps of the maximum of 
vorticity between the primary and the secondary cyclones. Such scenarios can 
shed light on the observed disruptions of the TC trajectories. 

On the basis of the present study we can conclude that the mcRSW model 
is able to capture some of essential dynamical features of TC and allows 
for semi-quantitative predictions of trajectories and intensities of TC at low 
computational cost. In particular, our investigation of interactions of TC- like 
vortices with topography reveals complex dynamical patterns and high sensitivity 
of the observed scenarios to various parameters, although the number of these 
latter is massively reduced as compared to “big” models. This shows, on the 
one hand, an intrinsic difficulty and complexity of the landfall forecasts but, on 
the other hand, a utility of simple models, because while a thorough analysis 
of landfall scenarios with mcRSW is affordable, this is hardly the case of “big” 
models. 

It is important to emphasize that the model can be further improved by 
including precipitable water in a liquid phase, together with vaporization, and its 
baroclinic two-layer version is also available Rostami and Zeitlin (2018). Radiative 
cooling and surface drag can be straightforwardly incorporated, as well, ensuring 
more realism. However, already the present results, obtained in the most simple 
version, clearly demonstrate the potential of the model. 

Let us finally stress that, independently of tentative utility for TC- related 
applications of the results in the moist-convective environment, their “dry” 
counterparts presented above contribute to better understanding of the dynamics 
of intense vortices, by providing new, and sometimes unexpected, information on 
their behaviour on the beta-plane and their interactions with topography, like, for 
example north-eastern deflection of trajectories, and swaps of vorticity maxima 
between primary and secondary vortices over the island. 
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